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LSD: Lyman-break galaxies Stellar populations and 
Dynamics. I: Mass, metallicity and gas at ^^3.1 * 
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ABSTRACT 

We present the first results of a project, LSD, aimed at obtaining spatially-resolved, 
near-infrared spectroscopy of a complete sample of Lyman-Break Galaxies at z~3. 
Deep observations with adaptive optics resulted in the detection of the main optical 
lines, such as [OIIJA3727, H/?, and [OIIIJA5007, which are used to study sizes, SFRs, 
morphologies, gas-phase metallicities, gas fractions and effective yields. Optical, near- 
IR and Spitzer/IRAC photometry is used to measure stellar mass. We obtain that 
morphologies are usually complex, with the presence of several peaks of emissions 
and companions that are not detected in broad-band images. Typical metallicities 
are 10-50% solar, with a strong evolution of the mass-metallicity relation from lower 
redshifts. Stellar masses, gas fraction, and evolutionary stages vary significantly among 
the galaxies, with less massive galaxies showing larger fractions of gas. In contrast with 
observations in the local universe, effective yields decrease with stellar mass and reach 
solar values at the low-mass end of the sample. This effect can be reproduced by gas 
infall with rates of the order of the SFRs. Outflows are present but are not needed 
to explain the mass-metallicity relation. We conclude that a large fraction of these 
galaxies are actively creating stars after major episodes of gas infall or merging. 

Key virords: Galaxies: abundances; Galaxies: formation; Galaxies: high- 
redshift; Galaxies: starburst 



1 INTRODUCTION 

Metallicity is one the most important properties of galaxies, 
and its study is able to shed light on the details of galaxy 
formation. It is an integrated property, not related to the 
present level of star formation, but rather to the whole past 
history of the galaxy. In particular, metallicity is sensitive to 
the fraction of baryonic mass already converted into stars, 
i.e., to the evolutionary stage of the galaxy. Also, metallicity 
is affected by the presence of infalls and outflows, i.e., by 



* Based on observations collected with ESO/VLT (proposals 
075.A-0300 and 076.A-0711), with the Italian TNG, operated 
by FGG (INAF) at the Spanish Observatorio del Roque de los 
Muchachos, and with the Spitzer Space Telescope, operated by 
JPL (Caltech) under a contract with NASA, 
t E-mail:filippo@arcetri. astro. it 



feedback processes and by the interplay between the forming 
galaxy and the intergalactic medium. 

The existence of a clear correlation between lumi- 
nosity and metallicity o f gala xies is known from the 
late '70s (jLegueux et all 1 19791 ). It is now clear that 
local galaxies follow a well-defined mass-metallicity rela- 
tion, where galaxies with larger stellar mass o r larger 
circular velocity have h i gher metallicities ([Garnetd 



iPerez-Gonzalez et al. I I2OO3I |Pilvugin_et al.l 12004: 



2008 



Tremonti et al.l |2004|; iLee et al] l2006l: iPanter et all 



2008 



.Kewlev fc Ellison _ _ . 

Michel-Dansa c et all l2'o08l : iLiu et al.l 120081 ). The 

of the relation is uncertain because several 



a iHavashi et all I2OO8I : 



origm 
effects can 

be, and probably are, active. It is well known that, in 
the local universe, starburst gala^xies eject a significant 
amount of metal-enriched gas into the intergalactic medium 
because of energetic feedback from exploding SNe (see, for 
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example, Lehnert fc HeckmarJ Il996al : I Mori et alj |2004| and 
IScannapieco et aL 20081 . and references therein). Numerous 
studies have s hown that outflows are a l ready present at 
high redshift s JPettini et al.1 1200 ll . |2002| : iFrve et all |2002| : 
IWeiner et al.l |2008| ). Outflows are expected to be more 
important in low-mass galaxies, where the gravitational 
potential is lower and a smaller fraction of gas is retained. 
As a consequence, higher mass galax ies are expecte d to 
be more meta l rich (see, for example, lEdmundill990l and 
lGarnett|[2003 ). 

A second possibile effect shaping the mass-metallicity 
relati on is related to the well kn o w effect of "downs izing" 
(e.g., iGavazzi fc Scodeggiol Il996l : ICowie et al.1 Il996l ). i.e., 
lower-mass galaxies form their stars later and on longer 
time scales than more massive systems. At a given time, 
lower mass galaxies have formed a smaller fraction of their 
stars, therefore ar e expected to show lower me tallicities 
l|EUison et al.ll2008l ). Recently. [Calura et al.l (|2009l ) have ex- 
plained the evolution of the mass-metallicity relation up to 
2=3.5 as due to an increase of the efficiency of star forma- 
tion with galaxy mass, without invoking differential galactic 
outflows. 

Other possibilities exist, for example some properties 
of star formation, as the initial mass functio n (IMF), could 
chan ge systematically with galaxy mass ijKoppen et al.l 
l2007h . 

All these effects have a deep impact on galaxy evolution, 
and the knowledge of their relative contributions is of cru- 
cial importance. Different models have been built to repro- 
duce the shape of the mass-metallicity relation in the local 
universe, and different assumptions produce divergent pre- 
dictions at high redshifts (z >2). It is possible to test these 
predictio ns as metallicities can now be measured at hi^ 
redshifts llKobuliiickv fc Koo 2000; Pcttini et al. 2001, 2003; 
iMaier et al.ll2006l ). To explore this issue several groups have 
observed the mass-m etallicity re lation in the distant uni- 
verse , around «=0.7 (ISavaglio et al. 20051: Rodri^ues et al.l 



l2008l ). betwee n and 1.5 j Cowic fc Bargcr' '2008') between 
0.5 a nd 1.2 llLamareille et aL 2008; Pcrcz-Montero et al.l 
20081). and at z~2 (|Erb et al.|[2006ai : iHaflidav et al.l l2008l : 



Havashi et al.|[2008l ) . They have found a clear evolution with 



cosmic time, with metallicity decreasing with in creasing red- 
shift, for a given stellar mass. lErb et al.l (|2006a 'b') used near- 
IR long-slit spectra of a large number of LBGs at z^2 to 
estimate the average amount of extinction and estimate the 
intrinsic star formation rate (SFR). By assuming that the 
Schmidt law holds at z=2 (see below), they also derive the 
fraction of baryonic mass in gas, and measure the effective 
yields. 

There are several reasons to explore even higher red- 
shifts. The redshift range at z^3-4 is particularly interest- 
ing: it is before th e peak of the cos mic star formation density 
(see, for example, iHopkins fc Bea com 2006; Mannucci ct a,l.! 
[2003), only a smaU fraction ('^15%. Pozzetti et al..,2007. ') of 
the total stellar mass (Af,) has already been created, is the 
redshift range where the most massiv e early- type galaxie s 
are expected to form (see, for example. ISaracco et al.ll2003l ). 
the number of merger s among the ga l axies i s much larger 
than at later times (jConselice et al.l l2007l : IStewart et al.l 
l2008l 'l. As a consequence, it is above z~3 that predictions of 
different models tend to diverge significantly, and it is im- 



portant to sample this redshift range observationally. The 
observations are really challenging because of the faintness 
of the targets and the precision required to obtain a reliable 
metallicity. Nevertheless, the new integral-field unit (IFU) 
instruments on 8-m class telescopes are sensitive enough to 
allow for the project. 

Metallicity at z~3 can be obtained by measuring 
the fluxes of the main optical emission lines ([OII]A3727, 
H/3, [OIII]A5007, Ha), whose ratios have be en calibrated 
again s t metallicity in the local universe (|Nagao et al.l 
I2OO6I : iKewlev fc EUisonI I2OO8I and references therein). Of 
course, this method can be applied only to line-emitting 
galaxies, i.e., to low-extinction, star-forming galaxies, 
whose lines c an be observed even at h i gh-red s hifts (e.g., 
Teplitz et all I2OO0I : iPettini et all I2OOII . |2002| : lErb et all 



2006al ). In contrast, the gas metallicity of more quiescent 



and /or dust extincted galax i es, lik e Extremely Red Galax- 
ies (ER Os, Mannucci et al.l I2002I '). Distant Red Galaxies 
(DRGs. iFranx et al.ll2003l ). and Sub-mm Galaxies (SMGs, 
" iOOSL 



IChapman et all 120051 ). cannot be easily measured at high 
redshifts. Stellar metalli cities, revealed by absorptio n lines, 
can also be measured ([Gallazzi et al.| [2006; Pantc i~et al.l 
' 2OO8I ; iHallidav et al.l 120081 ) . but at z >3 this is even more 
observationally challenging. 

Together with metallicity, the dynamical properties of 
galaxies have a special role for their understanding. Dynam- 
ics is directly related to the models of galaxy formation and 
it is the most direct way to probe the content of dark matter. 

Not much is known about dynamics of high red- 
shift galaxies as high-resolution, h igh-sensitivity, s p ectro - 
scop ic observati o ns ar e required. iNesvadba et al.l (|2006l ) 
and I St ark et al] (|2008l ) observed highly lensed LBGs at 
z~3.1 to obtain the velocity fleld. In these cases, the 
high spatial resolution on the source plane is provided 
by the presence of the gravitational lens, distorting the 
galaxy and stretching i t s app arent dimensions to several 
arcsecs. iNesvadba etal] l|2006D found rotation velocities of 
75 km/sec within the central lens-corrected 0.5 kpc, re- 
vealing the presen ce of a dynamical mass log(M/MQ)=9.3. 
IStark et al.l (|2008|) found a well ordered so urce, similar to lo- 
cal spirals. iForster Schreiber et al.l l|2006l) obtained seeing- 
limited near-IR integral field spectroscopy of a large sam- 
ple galaxies at 2~2. A significant fraction of this com- 
posite sample of spatially extended, Ha-emitting galaxies 
shows dynamical properties consistent with being the pre- 
cursor of normal disk galaxies i n the l ocal universe . Other 
rece nt works by iBouche et al.l (|2007l ). iGenzel et all (I2OO8I ) 
and ' Cresci et al.1 (|200^ are addressing this problem, and 
the amount of accumulated information is rapidly growing. 

Aimed at studying the high-redshift metallicity and dy- 
namics, we started two projects, AMAZE (Assessing the 
Mass- Abundance redshift Evolution) and LSD (Lyman- 
break galaxies Stellar populations and Dynamics j3. For both 
projects we observed a sample of galaxies at z—S-i with an 
IFU 3D spectrograph in order to derive their chemical and 
dynamical properties. 3D spectroscopy is a key technique 



^ This acronym is not to be confused with two others often used 
by astronomers, one o f whic h is "Lens Structure and Dynamics" 
bv lKoopmans fc Treul 1I2OO3I') . 
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Figure 1. Distribution in R magnitude and (G-R) color of 
the LSD targets (solid blue dots) compared to the spectroscop- 
ically confirmed targets of the Steidel et al. (2003) sample with 
2.5<z<3.3 (small dots). The solid blue histograms show the distri- 
butions of these quantities for the LSD sample, which are similar 
to those for the total parent sample (dashed histograms). 

for these studies as it allows to derive the full velocity field 
of the galaxies, without the need to restrict the study to a 
slit. Also, integrated spectra of the galaxies can be obtained 
without slit losses. This is particularly important when the 
flux ratios between lines at different wavelengths are derived, 
as differential losses can spoil the results. 

For AMAZE, described in iMaiolino et all (|2008l ). we se- 
lected 30 galaxies with deep Spitzer/IRAC photometry (3.6- 
8 fj,m), an important piece of information to derive a reliable 
stellar mass. These galaxies were observed, in seeing-limited 
mode, with the IFU spectrometer SlNFONl on ESO/VLT. 

The aim of LSD is not only to measure metallicities, but 
also to obtain spatially-resolved spectra to measure dynam- 
ics and spectral gradients. For this reason we used adaptive 
optics to improve the spatial resolution. Also, we selected 
a complete, albeit small, sample of LBGs, in order to ob- 
tain results that can be applied to the LBG population as 
a whole. This is the first paper about LSD, focusing on the 
chemical properties, while the dynamical properties will be 
examined in a forthcoming paper. 

In the following, we adopt a ACDM cosmology with 
Ho=70 km/sec, flm~0.3 and 57a=0.7. At z=3.1, the uni- 
verse is 2 Gyr old, 15% of its current age. 



2 SAMPLE SELECTION 

Our s ample of LBGs was extracted from the ISteidel et al.l 
catalog, which contains about 1000 spectroscopically- 
confirmed LBGs, selected in 17 fields down to a limiting 
magnitude of Rab ~ 25.5. The parent catalog was searched 
for the presence of bright stars (R<15) within ~35" of each 



source, resulting in the selection of ~100 sources. The pres- 
ence of a bright foreground star is needed to drive the adap- 
tive optics module and obtain a resolution comparable with 
the difi'raction limit of a 8m class telescope, about 0.1" in 
the near-IR. The best combinations of star brightness and 
distance from the target were selected to obtain the final 
sample of 10 objects presented in Table [1] 

We emphasize that no other property was used to select 
the targets, thus the present sample, albeit sm all, should be 
repres entative of the LBG population in the ISteidel et al.l 
(2003) sample. This is shown in figure [1] which compares 
the distributions of R-band magnitude and (G-R) color of 
the selected targets with those of the parent population of 
LBGs. The distributions of both quantities are similar, with 
no apparent biases. 

In order to check the possible presence of Active Galac- 
tic Nuclei, we inspected the available X-ray observations 
of these galaxies. For two objects, DSF2237b-D28 and 
DSF2237b-MD19, no X-ray observation is available. The 
other sources lay in the field of a deep Chandra observa- 
tion (SSA22) and in that of a quasar observed with XMM- 
Newton. No source is detected in the 0.5-10 keV range. In all 
cases, the off-center position of the sources, and/or the short 
exposure time imply a relatively high upper limit for the X- 
ray fiux, which, given the high redshift of the sources, cor- 
responds to X-ray luminosities larger than ~ 10*"^ erg s~^. 
Although these limits cannot rule out the presence of an 
AGN, there is no direct evidence for it. In all cases we can 
exclude a dominant, quasar-like AGN contribution. 



3 OBSERVATIONS AND DATA REDUCTION 

The selected objects have been observed with several instru- 
ments in order to obtain a complete data set. 

3.1 Integral field spectroscopy 

The targets were ob served at ESO/VLT by using the 
SINF ONI instrument (lEisenhauer e t al. 20o3; teonnet et al] 
|2004 ). an IFU spectrograph fed by an adaptive optics mod- 
ule with a Natural Guide Star (NGS). The IFU splits the 
incoming light into several adjacent slits, and a spectrum 
is obtained for each location in the field-of-view. Several 
spatial scales are available, corresponding to different spa- 
tial resolutions. For 8 of the 10 targets we used the 50x100 
mas/pixel scale, providing a field-of-view of 3"x3". Such a 
scale is a good match of the theoretical diffraction limit in 
the K-band on the VLT telescope, which is ~70 mas, and the 
field-of-view is appropriate to observe galaxies of dimensions 
of ~1". Two objects constitute a pair of interacting galaxies 
at a distance of about 2", in this case we used a larger scale, 
125x250 mas/pixels with a total field-of-view of 8"x8". 

In all cases we used the H-l-K grating, providing a simul- 
taneous coverage of these two near-IR bands (1.5-2.4/^m), 
with resolution A/AA ~1400 in the H band and ~1800 in 
the K band. 

Observations were divided into independent sequences 
of about 45 min each, due to the Hmitation of the ESO 
service observing strategy based on observing blocks (OBs) 
of about 1 hour of total time. During each OB, an ABBA 
nodding was used, putting the target in two positions of the 
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Table 1. Target list and summary of the observations 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Object 


R.A. (J2000) Dec. 


Rab 






SINFONI 






SPITZER 














Dist 


R mag 


Scale 


Texp 


T(l+3) 


T{2+4) 


N/A 


SSA22a-C30 


22:17:19.3 


+00:15:44.7 


24.2 


3.103 


35" 


13.0 


50x100 


240 


50 


50 


A 


SSA22a-C6 


22:17:40.9 


+00:11:26.0 


23.4 


3.097 


35" 


12.0 


125x250 


280 


83 


133 


N 


SSA22a-M4 


22:17:40.9 


+00:11:27.9 


24.8 


3.098 


35" 


12.0 


125x250 


280 


83 


133 


N 


SSA22b-C5 


22:17:47.1 


+00:04:25.7 


25.2 


3.112 


11" 


14.9 


50x100 


240 


83 


133 


N 


DSF2237b-D28 


22:39:20.2 


+11 55 11.3 


24.5 


2.933 


26" 


13.1 


50x100 


240 


96 


96 


A 


DSF2237b-MD19 


22:39:21.1 


+11 48 27.7 


24.5 


2.611 


20" 


14.1 


50x100 


200 








Q0302-C131 


03:04:35.0 


-00:11:18.3 


24.5 


3.235 


26" 


13.4 


50x100 


240 


332 


133 


N 


Q0302-M80 


03:04:45.7 


-00:13:40.6 


24.1 


3.414 


15" 


14.5 


50x100 


240 


166 


399 


N 


Q0302-C171 


03:04:44.3 


-00:08:23.2 


24.6 


3.328 


33" 


14.6 


50x100 


240 


332 




N 


Q0302-MD287 


03:04:52.8 


-00:09:54.6 


24.8 


2.395 


13" 


15.0 


50x100 


160 


249 




N 



Columns: (1-3): Target name and coordinates; (4): Target R-band magnitude; (5): Redshift from emission lines; (6-7): Guide star 
distance and R-band magnitude; (8-9): spatial scale (mas) and total exposure time (min) of the SINFONI observations; (10-11): 
exposure time (min) for Spitzer/IRAC observations in channels 1+3 and 2+4; (12) New (N) or Archive (A) Spitzer/IRAC data. 



field-of-view about 1.5" apart. This is a very efficient way 
to observe compact objects as no time is lost in observing 
the sky. Nevertheless, in principle this could introduce some 
self-subtraction of the wings of the targets more extended 
than 1.5". 

A maximum seeing of 0.8" was requested to carry out 
the observations, and an image of the reference star was ob- 
served at the beginning of each OB to monitor the conditions 
of the atmosphere during the observations. 

Data were reduc ed by using the ESO pipeline 
l|Modigliani et all |2007|) w i th th e improved sky subtrac- 
tion described by Davied |20o3). After flat-fielding, sky- 
subtraction, correction for distortions and wavelength cali- 
bration, the 2D data of each OB are mapped into a 3D data 
cube with dimensions 32 x 64 x 2048 pixels. The cubes of 
the different OBs are then combined together by measuring 
the relative ofi'sets from the detected centroid of emission. 
We note that using the coordinates in the image header of- 
ten does not provide the requested precision (<0.l"), as the 
measured pointing uncertainties are of the order of 0.2". 



3.2 Photometric and morphological data 

The orig inal dataset of opt ical photometry in the UGR 
bands bv lSteidel et al.l (|2003l ) is available for all the objects. 
Such a limited wavelength range does not allow to obtain 
reliable stellar masses by fitting the spectral energy dis- 
tribution (SED). The optical data must be complemented 
by near-IR photometry. J and K-band photo metry of 4 
objects was published by IShaplev et al] l|200ll ). We have 
observed the remaining 6 objects in the J and Ks filters 
with the NICS camera (|Baffa et al.l I2OOII ) on the Itahan 
3.6m telescope TNG. Exposure times of 1-1.5 hours in J 
and 2-3.5 hours in Ks were used, according to the expected 
near-IR magnitudes. NICS data were reduced by the 
automatic pipeline SNAF0. 

Photometry at longer wavelength s (3.6-8Atm) is avail- 
able in the Spitzer/IRAC archive l|Werner et all |2004 



http: / / www. arcetri . astro, it / ~filippo / snap 



iFazio et al.ll2004l ) for some of the objects in the SSA22 field. 
For the remaining objects we obtained targeted observation 
with Spitzer/IRAC during cycle 7, as listed in Table [1] We 
used the post-BSD products of the pipeline, and measured 
aperture photometry. 

Finally, HST will observe these objects after SM4, both 
in the optical with ACS and in the near-IR with WFC3, pro- 
viding accurate photometry and broad-band morphologies. 
In conclusion, these objects will have a complete set of pho- 
tometric and spectroscopic observations with detailed mor- 
phological information at optical and near-IR bands, and 
will be suitable for detailed studies. 



4 RESULTS AND MEASURED QUANTITIES 
4.1 The detected lines 

The final SINFONI data cubes were searched for the po- 
sition of the brightest lines, either [OIIIJA5007 or Ha, to 
define the spatial position of the target. The main optical 
lines ([OII]A3727, H/3, [OIII]A5007 and Ha) are detected in 9 
out of 10 galaxies, the only exception being Q0302-MD287. 
The reasons for this non detection are not clear. It is possi- 
ble that its emission lines fall below our detection threshold, 
or that they coincide with bright sky lines. The SED fitting 
described in sect. l4.3l points toward the existence of a large 
amunt of dust in this object, with Av between 0.9 and 1.1 
mag. As a consequence, it is also possible that the opti- 
cal lines are not detected because of the presence of large 
amount of dust extinction. 

The choice of the aperture for the extraction of the 
spectra is a important step. Measuring metallicity requires a 
good knowledge of line ratios, which are best obtained inside 
small apertures on the central regions of the galaxies, where 
the signal-to-noise is higher. In contrast, measuring total 
SFRs requires larger apertures to maximize the fraction of 
line emission recovered by aperture photometry. To measure 
all these quantities, we extracted complete spectra inside 
a circular aperture of 0.35" of diameter, corresponding to 
7 pixels, and an aperture correction for a point source was 
applied to recover the emitted fiux of the central part of 



LSD: Mass, metallicity and gas at z^S.l 5 



Table 2. Properties of the detected lines 



(1) 

Object 


(2) 
lino 


(3) 

z 


(4) 
Flux(") 


(5) 
FactC") 


SSA22a-C30 


[Oil] 3727 

H/3 
[OIII]4958 
[GUI] 5007 


3.1024 
3.1030 
3.1034 
3.1026 


0.30 ± 0.12 
0.26 ± 0.10 
0.34 ± 0.10 
1.28 ± 0.22 


3.18 


SSA22a-C6 


[Oil] 3727 

[OIII]4958 
[GUI] 5007 


3.0966 
3.0970 
3.0966 
3.0968 


1.21 ± 0.38 
0.76 ± 0.45 
2.26 ± 0.78 
5.50 ± 0.62 


1.09 


SSA22a-M4 


[Oil] 3727 

[GIII]4958 
[GUI] 5007 


3.0961 
3.0976 
3.0965 
3.0978 


1.46 ± 0.23 
0.69 ± 0.25 
1.08 ± 0.30 
3.66 ± 0.49 


1.06 


SSA22b-C5 


[NoIII]3869 

[GIII]4958 
[GUI] 5007 


3.1098 
3.1122 
3.1120 
3.1121 


0.38 ± 0.14 
0.50 ± 0.14 
1.01 ± 0.22 
3.28 ± 0.24 


1.57 


D2237b-D28 


[Oil] 3727 
[GIII]4958 
[GUI] 5007 


2.9352 
2.9328 
2.9324 


0.74 ± 0.16 
0.50 ± 0.16 
1.86 ± 0.36 


1.84 


D2237b-MD19 


Ha 


2.6093 
2.6108 


0.60 ± 0.16 
2.96 ± 0.46 


1.95 


Q0302-C131 


[Oil] 3727 

tip 

[GIII]4958 
[GUI] 5007 


3.2349 

3.2341 
3.2346 


0.80 ± 0.32 
U.4D it U.iD 
0.90 ± 0.20 
2.62 ± 0.52 


1.29 


Q0302-M80 


[OII]3727 

H/3 
[GIII]4958 
[GUI] 5007 


3.4136 
3.4138 
3.3990 
3.4137 


0.56 ± 0.18 
0.40 ± 0.12 
0.42 ± 0.20 
1.32 ± 0.26 


1.23 


Q0302-C171 


[OII]A3727 

H/3 
[OIII]A4958 
[OIII]A5007 


3.3287 
3.3289 
3.3280 
3.3279 


0.56 ± 0.24 
0.18 ± 0.08 
0.28 ± 0.08 
1.10 ± 0.12 


1.92 



Columns. (4): line flux inside an aperture of 0.35" with correction 



for the PSF, 
total flux 



units of 10 erg cm sec ; (5): scaling factor to 



the galaxies. We chose a circular aperture, instead of more 
complex apertures based on object morphology and surface 
brightness, in order to avoid any problem with possible 
variations of the PSF with wavelength and a non perfect 
spatial alignment of the spectra. The fraction of total flux 
recovered by this procedure was estimated by computing 
the curve of grow of the photometry of the brightest line as 
a function of aperture radius. A correction for the missing 
flux is applied when the "total" quantities are needed. The 
results are shown in Table (5] listing the line fluxes inside the 
circular aperture and the correcting factor to total fluxes. 
Except for one galaxy, SSA22a-C30, composed by several 
clumps of emission, for all the other galaxies the tabulated 
fluxes comprise more then half of the total flux. 



4.2 Morphologies and sizes 

Fig. [2] shows the images of all the targets. Most of the 
galaxies show a sharp peak of line emission and some 
extended structure. In several cases, the objects appear to 
be composed by several clumps of emission. The distance 
between clumps is usually comparable to the dimension 
of the objects as revealed by ground-based imaging. In a 
few cases the distribution of the secondary peaks follow 
the continuum light, while in some cases we detected 
faint "companions" that are not detected in broad-band. 
VLT has a small pointing error, of the order of 0.2", but, 
compared to our resolution, this is enough to introduce a 
signiflcant uncertainty in the relative astrometry between 
UV and line images. As a consequence, we cannot be 
sure where the line emission originates with respect to 
the UV image. Assuming that the main peak of line 
emission corresponds to the maximum of the UV light, 
the line-emission "companions" could be both secondary 
peak of star formation and line emission inside a single 
formed galaxy, or the signature of an ongoing merging. 
A detailed comparison with broad-band, rest-frame-UV 
morphology will be given on the basis of future HST images. 

Images obtained wit h a NGS usually ha ve a non-circular 
PSF (see, for example. ICresci et al.|[2005h . Increasing the 
distance of the scientiflc target from the NGS, the PSF be- 
comes larger. The variation in the radial direction, i.e., to- 
wards the NGS, is larger than in the tangential direction. 
As a result, when this effect is dominant, the objects appear 
to be elongated toward the NGS. The magnitude of the ef- 
fect depends on the conditions of the atmosphere and on the 
distance of the target from the NGS, and is of the order of 
0.l"-0.3". In our images, elongation toward the NGS is seen 
only in one case, meaning that in most cases the atmospheric 
conditions were good enough for AO observations and the 
isoplanatic angle in the K band was not signiflcantly smaller 
than the separation between the target and the guide star. 

No bright point sources are present in our field-of-view, 
preventing us from measuring the flnal PSF on the com- 
bined data themselves. To estimate the PSF we used the 
on-axis image of the NGS taken every hour at the beginning 
of each OB. The on-axis PSF FWHM was degraded to take 
into account the distance of the target galaxy from the 
NGS. For this we used the recipes in Cresci et al, (.20051 ). 
using the values of the isoplanatic angle provided by the 
ESO database, and averaging over the different values 
obtained for the different OBs. As the seeing was always 
comprised in the narrow range 0.6-1.0" FWHM, no large 
differences are seen in most of the OBs. In all cases we 
obtained values of the PSF FWHM of the order of 0.2"for 
the fields observed with the smallest spatial scale, and 
0.3"for the fields with the largest scale. These values are in 
good agreement with the results of the simulations provided 
by ESO, and were used in the following analysis. 

Galaxy sizes have been estimated by fitting a Moffat 
function to the image in the brightest line, and the results 
have been deconvolved for the PSF. The measured half- 
light radii are listed in Table |31 Large uncertainties on the 
resulting sizes are present, both because images usually 
have low signal-to-noise ratio (SNR) and because the PSF 
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Figure 2. Images and spectra of the targets. For each object we show , from lef t to ri ght: morphology in the emission line (Ha for 
MD19, [OIII]A5007 for the others), ground-based R-band image by Stei del et all l|2003h . and galaxy spectrum compared to the sky 
spectrum. The horizontal bars in the line images are 1" long. First column, from top to bottom: SSA22a-C30, the interacting system 
composed by SSA22a-C6 (south, upper spectrum) and SSA22a-M4 (north, lower spectrum), SSA22b-C5, DSF2237b-D28. Second column: 
DSF2237b-MD19, Q0302-C131, Q0302-M80 and Q0302-C171. 



is not accurately known. Also, companions and secondary 
peaks of emission are not considered in measuring the 
dimensions. As this is also true for line fluxes, the derived 
metallicities (section I4.6|l and gas mass (section I6.3|l . are 
consistently derived from the main peak only. 

Since we are observing faint galaxies with high-spatial 



resolution, surface brightness is probably our limiting fac- 
tor. The dimensions of the detected line-emitting regions 
are sometimes significantly smaller than the dimensions ob- 
tained in the UV by ground-based imaging, even taking into 
account the larger PSF. They are also smaller than object 
sizes obtained with SINFONI/VLT without the use of AO by 
the Spectroscopic Imaging survey in the Near-infrared with 



LSD: Mass, metallicity and gas at z^S.l 7 



SINFONI (SINS) jForster Schreiber et alfeood . although at 
zr^ 2.2) and AMAZE (|Maiolino et al.ll2008h collaborations. 
We interpret these results as evidence that our detected sizes 
are indeed limited by surface brightness. The average lim- 
iting surface brightness is 1.2 x 10~^® erg/cm'^ /sec/arcsec^. 
Such a sensitivity can be translated into a limiting SFR den- 
sity at 2:~3.1 detected with an emission line. For H/3, the 
limiting SFR density is ~4/M0/yr/kpc^; for [OIII]A5007, 
assuming an average ratio of [OIII]A5007/H/3=4, the mini- 
mum detectable SFR density is ^I Ma/yr/kpc^. Thi s is sim- 
ilar to the SFR density at z^2.2 in lErb et all l|2006d ). and it 
is adeq uate to detect a significant fraction of the local star- 
bursts l|Kennicuttll998l ) . Nevertheless, the outer, low surface 
brightness parts of these galaxies are probably lost in the 
noise, and morphologies could appear more compact than 
they actually are. 

The measured half-ligh radii range between 0.7 and 
2.4 kpc, with a median value of 1.36 kpc. These radii 
ar similar to those in the s ample of local starbursts 
in iLehnert fc HeckmanI l|l996bl ). who measure a median 
value of 1.70 kpc. Considering that our observations are 
less sensitive to the oute r parts of the galaxies than 
iLehnert fc HeckmanI (Il996al ). we conclude that LBGs tend 
to be as extended as local starbursts. 



4.3 The Stellar mass 

Stellar mass es (M<.) are derived by fitting the SED with 
Hyperzmass jPozzetti et akllioOTl ). a modified version of th e 
photometric redshift code Hyperz iBolzoncUa ct al1 l2000h . 
These fits also provide estimates of age, current SFR and 
dust extinction of the dominant stellar population. The 
presence of good IRAC photometry allows the determina- 
tion of reliable M, as the SED is sampled up to the rest - 
frame J band. We have used the lBruzual fc CharlotI (|2003l ) 
spectrophotometric models of galaxy evolution and smooth 
exponentially decreasing SFHs, constraining the age to be 
s maller than the Hubbl e time at galaxy redshift, as detailed 
in lPozzetti et all (|2007l ). 

The photometric stellar masses have typical dispersions 
due to statistical uncertainties and degeneracies of the or- 
der of 0.2 dex. Metallicity has little influence on the SED, 
with uncertainties of ~0.05 dex, and we have set it to Z=0.2 
Zq, in agreement with what derived from the emission lines 
as described below. The addition of secondary bursts to a 
continuous star formation history produces systematically 
higher (up to 40% on average) stellar masses, while popula- 
tion synthesis models with TP-AGB phase (Maraston 2005) 
produce a systematic shift of ~0.1 dex towards lower M«. 
Finally, the uncertainty on the absolute value of the M« due 
to the assumptions on the IMF is within a factor of 2 for the 
typ ical IMFs usually adopted in the literature. We have used 
the lChabrierl (2003) IMF with lower and upper cutoffs of 0.1 
and 100 Mq. The resulting stellar masses can be scaled to 
the sta ndard Salpeter IMF by multiplying them by a factor 
of 1.7 (jPozzetti et al.ll2007l ). 

The results are shown in Table O The listed errors 
are due to both statistical uncertainties and degeneracies. 
Masses, expressed in log(M»/M0), range from 8.96 to 10.33, 
with a log averag e of 9.82. Thes e value s are consistent with 
those derived by IShaplev et al.l (|200ll ) for their sample of 
LBGs at z~3, averaging at 10.05 when converting their val- 



ues to our IMF. Lower redsh ift BM/BX galaxies in the z~2 
sample of lErb et all (|2006bh are also quite similar, with an 
average log(M*/M0) of 10.22, within ~0.3 of our average 
value. 

Galaxy age, current SFR and dust extinction on the 
continuum Av,c cannot be individually well constrained be- 
cause they suffer of significant degeneracy. The errors asso- 
ciated to each of these three quantities are large, and this 
is particularly important for extinction which is used to es- 
timate the total SFR from emission lines. The range of the 
allowed values of Av,c is shown in Table[3l Despite the large 
errors, a possible correlation is seen between stellar mass and 
age, i.e., less massive galaxies are best-fitted by younger tem- 
plates. The best-fitting age is of the order of 10^'^ yr for the 
less massive galaxies, and about 10^ yr for the most massive 



4.4 Photoionization conditions 

The values of SFR and metallicity derived in the next sec- 
tions are correct only if the optical line emission is domi- 
nated by star formation and the conditions are not too dif- 
ferent from those in the local universe. Even if the observed 
line ratios are typical of starburst galaxies, several problems 
could be present. 

The first concern is that the presence of an AGN could 
alter line ratios and produce a spurious value of the inferred 
metallicity. No evidence of AGN is seen in the emission line 
shape, where the width of the forbidden lines is consistent 
with those of the permitted lines. This excludes the pres- 
ence of a dominant broad-line AGN, but it does not help in 
excluding a narrow- line AGN (see the discussion in sec. [5}. 

The second concern is that significant differences in the 
conditions of star forming regions can be present, as LBGs 
at z^3 have larger SFRs than most galaxies in the local 
universe. In particular, any evolution of the photoioniza- 
tion conditions could mimic a change of metallicity. The 
presence of evolution can be studied by comparing the 
flux ratios [NII]/Ha vs. [OIII]/H/3, and several studies at 
lower redshift indicate that such an evol u tion actually exists 
(Shaplcv ct al. 20qi lErb et all l2006al : iBrinchmann et all 
I^OOSl : iLiu et al.ll200a ). Line ratios discrepant from the lo- 
cal relation are often associated with higher SFR surface 
densities, interstellar pressures, and ionization parameters. 

We have verifled that all the observed line ratios are 
fully consistent with excitation from hot stars, but we 
cannot study the presence of evolution because in most 
cases we cannot observe Hq. Neverthele ss, several authors 
(|Brinchmann et al.ll2008l : iLiu et al.ll2008l ) conclude that the 
influence of different ionization conditions on the derived 
metallicities is likely to be low, of the order of 0.1 dex. 



4.5 SFR from optical lines and SED fitting 

SFR can be derived from the Ha and H/3 lines. Several 
conversion are possible, depending on the available data 
and physical condition of the galaxies. Similar to what 
is often done by many aut hors (e.g., lErb et al] l2006al : 
iForster Schreiber et al] I2OO9I) . we have used the "classi- 



cal" conversion factor by KennicuttI (Il998|). sc aling down 
the results by a factor 1.7 i Pozzetti et al. 20071 ) to convert 
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Table 3. Properties of the target galaxies 



(1) 

Object 



(2) 

log(M*/MQ) 



(3) 
range 



(4) 

SFRem 

(Me/yr) 



(5) 
(kpc) 



(6) 

12+lg(0/H) 



(7) 

lg{M9„,/M0) 



(8) 



(9) 

ig(ye//) 



SSA22a-C30 

SSA22a-C6 

SSA22a-M4 

SSA22b-C5 

DSF22-D28 

DSF22-MD19 

Q0302-C131 

Q0302-M80 

Q0302-C171 



1 n 

±U.OO_g 3g 

,+0.34 
-0.13 



9.4i: 



8.96 
9.78 



+0.38 
-0.22 
+0.28 
0.29 

lo.oetoir 

10 09+"'^° 
1 n n7+'''2^ 

lU.U/_g j^g 

10.06t«-° 



0.0-0.8 
0.3-0.5 
0.0-0.8 
0.0-0.5 
0.0-0.6 
0.0-1.2 
0.0-0.3 
0.0-0.6 
0.0-0.2 



29 
23" 



+81 
-21 
-11 



20l« 

18 
8 
+ 77 
-26 
+6 
-4 
17 



14' 
39 

10 
13 

5 



+ 2 
2 



1.48±0.44 
1.75±0.22 
2.01±0.27 
1.28±0.38 
1.50±0.48 
2.00±0.55 
1.27±0.37 
0.75±0.24 
1.25±0.39 



8.16 
7.95 
8.12 
7. 

8.20 

8.00 
8.36 
8.14' 



+0.20 
0.60 
+0.20 
0.50 
+0.25 
0.45 
+0.20 
-0.20 
+0.10 
0.25 



+0.25 
0.40 
+0.15 
0.15 
+0.25 
0.45 



9.66±0.41 
9.96±0.13 
9.96±0.34 
9.64±0.21 
9.62±0.25 
9.98±0.33 
9.57±0.14 
9.54±0.25 
9.26±0.11 



0.18 
0.66 
0.78 
0.83 
0.61 
0.69 
0.42 
0.36 
0.30 



+0.26 
0.12 
+0.08 
0.12 
+0.12 
0.28 
+0.08 
0.19 
+0.21 
0.20 
+0.25 
0.26 
+0.13 
0.07 
+0.16 
0.12 
+0.08 
0.04 



-2.63 



+0.2+0.3 
-0.7-0.2 
r, 99+0.2+0.1 
^•^^-0.4-0.2 
-1 84+0.2+0.4 
-^•°^-0.4-0.4 
7+0.2+0.3 
0.2-0.4 
+0.1+0.3 



-2.17"' 
-2.32 



_9 70 + 0.2+0.2 
'^•'''-0.4-0.1 
_9 07+0.1+0.2 
^•■^'-0.1-0.2 
1+0.2+0.1 
0.4-0.1 



-2.71"* 



Columns. (2-3) Stellar mass and dust extinction from SED fitting; (4): total SFR from emission lines; (5): half-light radius of 
line emission; (6): gas-phase metallicity (7-8): gas mass and gas fraction; (9): effective yields; the first error is due to metallicity, the 

second one to gas fraction 
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Figure 3. Comparison between the SFRs derived from SED 
fitting and line emission. Solid and empty dots show LSD and 
AMAZE galaxies, respectively. The dashed line shows equal 
SFRs. A good agreement is found between the two estimates of 
the SFR. 



them to the IChabrierl ()2003l ) IMF. The iKennicutd (|l998l ) 
conversion factor is based on the H/3 flux corrected for 
dust extinction, and therefore requires a good knowledge 
of the amo unt of extinction A v.e.i suffered by the emis- 
sion lines. ICalzetti et al.l (|2000l ) found that this is propor- 
tional to the extinction suffered by the continuum Av.c, 
as r neasured by SED fitting, wit h Av,e! = j4v,c/0.44 (see 
also iForster Schreiber et allboogl ). and here we have used 
this assumption. For DSF2237b-D28 the H/3 flux is derived 
from [OIIIJA5007 and the expected line ratio for the value 
of metallicity measured in the next section. The values of 
SFRs are listed in Table |3] and plotted in Fig. [3] 

As explained in the previous section, large ranges of 



extinction are allowed by our SED fitting, therefore the re- 
sulting SFRs have large errors, with typical uncertainties 
of one order of magnitude. This is common to most of the 
studies of high-redshift star-forming galaxies, even if some- 
times the errors on the SFR are not accurately discussed. 
Unfortunately, other possible conversion factors that could 
provide smaller uncertainties cannot be applied to high- 
redshift starburt galaxies. For example, several other esti- 
mates have been propo sed b viArgence fc LamareillQ (|200^ , 
iMoustakas et al.l (|2006| ) and IWeiner et al.1 l|2007h . based on 
a combination of two emission lines (such as [OII]A3727and 
H/3), on lines and broad-band photometry, or line fiux with 
no correction for dust extinction. These recipes are based on 
the empirical correlation among mass, SFR, dust extinction, 
luminosity and metallicity observed in local galaxies, and it 
is not possible to apply them for galaxies with very different 
properties. 

The SFR derived from the emission lines can be com- 
pared with the value resulting from the fit of the SED (see 
Fig-E}. Large uncertainties are present due to the errors on 
H/3 fiux, the uncertainties in dust extinction correction for 
both estimates of SFR, the intrinsic spread of both calibra- 
tions. Nevertheless a good agreement is seen between the 
two estimators of the SFR. 



4.6 Metallicity 

Gas-phase metallicities were derived by comparing the ob - 
serv ed line ratios wi t h the calibrations in lNagao et al.l (|2006l ) 
and lMaiolino et ahl (|2008l ). Metallicities are derived by a si- 
m ultaneous fit of all th e available line ratios, as explained 
in lMaiolino et af] (|2008l ) . Both the uncertainties in the ob- 
served line ratios and the spread of the calibration were 
considered, and in some cases the latter is the dominant 
contribution to the metallicity error. In practice, the re- 
sult is dominated by the R23 indicator or, similarly, by the 
[OIII]A5007/H/3 ratio, while the [OIII]A5007/[OII]A3727 ra- 
tio is used to discriminate between the two possible branches 
at low- and high-metaUicity. The [Nelll] 3869/ [Oil] 3727 fine 
ratio is also a sensitive indicator, when both lines are de- 
tected. The uncertainties on metallicity are of the order 
of 0.2-0.3 dex, and tend to be larger when the metal- 
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licity has values 12+log(0/H)~8, when both R23 and 
[OIII]A4958,5007/H/3 have a maximum. In this case a smaU 
uncertainty in the hne ratio produces a large uncertainty in 
metallicity. The results are listed in Table [S] LSD galaxies 
show metallic ities about 10%-50% solar, similar to the val- 
ues found by iPettini et al.l l|200ll ) for a different sample of 
LBGs at z~3. 

Dust extinction is not strongly affecting the re- 
sults as the only extinction-sensitive line ratio is 
fOIII]A5007 /[OII]A3727. This eff ect can be clearly seen in 
Figure 6 of iMaiolino et af] (|2008l ). showing that metallicity 
is barely dependent on extinction and, as a consequence, the 
observed line ratios leave extinction basically unconstrained. 
For this reason we have chosen to limit the extinction to be 
within the range allowed by the SED fitting. Leaving extinc- 
tion totally unconstrained does not affect the best-fitting 
value of metallicity but only the range of confidence. 



5 NOTES ON INDIVIDUAL OBJECTS 
5.1 SSA22a-C30 

This is the galaxy with the largest stellar mass. While H/3 
and [OIII]A5007 are in clean parts of the atmospheric spec- 
trum, [OII]A3727 and [OIII]A4958 are affected by OH lines. 
[OIII]A5007 emission is characterized by the presence of a 
secondary peak, about 1" apart from the main one. Other 
lower-luminosity peaks can be present in between, although 
with lower significance. As a result, the object has a com- 
plex, very disturbed morphology. In the ground-based UV 
image, this galaxy appears less disturbed. It is elongated as 
the line-emitting blobs and can comprise them all. 



5.2 SSA22a-C6 and SSA22a-M4 

This is a pair of interacting galaxies at a projected distance 
of about 2" and with velocity difference of about 300 km/sec. 
All hues, except [NeIII]A3869, fall far from the OH atmo- 
spheric lines, allowing for a good measurement of the metal- 
licity. C6, the southernmost object, is the brightest mem- 
ber both in UV and in line emission. It has an elongated 
structure in the UV, and also the line image shows some 
secondary peaks coincident with the UV elongation. M4 is 
compact in both UV and line images. 



5.3 SSA22b-C5 

This is the object with the smallest stellar mass, the only 
one with log(M./Mo)<9. While [OIII]A4958,5007 are well 
detected, [OIIJA3727 is hidden by a sky line. A marginally 
significant structure is detected near the expected posi- 
tion of the [NeIII]A3869 line, providing a value of the 
redshift (2=3.116) slightly larger that the value from the 
[OIII]A5007 line (z=3.112). This could be a real detection 
of the [NeIII]A3869 line, but both the low statistical signifi- 
cance and the small velocity offset make it quite uncertain. 
Here we assume that this is not a real line and do not use it 
to derive metallicity. The galaxy shape is regular and almost 
circular, with no significant structure around the main core. 



5.4 DSF22b-D28 

The line image is characterized by a compact nucleus sur- 
rounded by an extended faint emission, with a size much 
more compact than in the broad-band image. Its emission 
lines are quite faint and H/3 is not detected. 

5.5 DSF22b-MD19 

This is the galaxy of our sample with the lowest redshift, 
2=2.616, and it is the only one where Hq can be observed. 
This line falls in the K band at 2.37/^m, where the thermal 
contribution to the background is significant. Nevertheless, 
it is very well detected. The nearby [Nil] lines are not de- 
tected and we can put a limit to their fiux ratio with Ha 
of [NH]6583/Ha<0.3. Such a value excludes that the ioniza- 
tion could be dominated by an AGN. The [OII]A3727 line is 
outside the observed range and the [OIIIJA5007 line is barely 
detected as it is coincident with a sky line. As a consequence 
we cannot measure the metallicity of this object. H/3 is well 
detected in a clean part of the spectrum, and the Hq/H/3 
ratio is 4.9±1.5, significantly above the case-B value of 2.85. 
Assuming that this difference is due to a screen of dust fol- 
lowing the ICardelli et al.l (|l989l ) extinction law, we derive 
A-v,ei = 1.7. This is consistent with what is derived by the 
spectral fitting, Av,c=i-l- MD19 is very extended both in 
UV and in Ha. The stellar continuum of this object is clearly 
detected and appears to be displaced from the peak of the 
lines. 



5.6 Q0302-C131 

This galaxy is elongated both in line and UV images, with a 
shape reminding a classical edge-on disk. The position angle 
of the elongation in the line image is fully consistent with 
that in the UV image. Most of the lines are in a clean part 
of the spectrum, only the [NeIII]A3869 line is hidden by a 
sky line. 

5.7 Q0302-M80 

In the broad-band image this galaxy is very extended, with 
an elongated structure. In contrast, the line image appears 
very compact and regular. Both [Oil] A3727 and [Nelll] A3869 
fall near bright sky lines and cannot be detected. 

5.8 Q0302-C171 

This is the object with the faintest [OHI]A5007. While most 
of the main emission lines are not strongly affected by sky 
lines, the flux of [OH]A3727 is very uncertain because the 
line falls near a bright sky line at 1.613/im. 

5.9 Composite spectrum 

The average spectral properties of the LSD galaxies can 
be better shown by constructing an average spectrum. 
Also, fainter spectral features, below detection in any sin- 
gle spectrum, could be revealed by summing up all the 
spectra. Fig. |4] shows the average, composite spectrum of 
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Figure 4. Composite spectrum obtained by averaging the spectra of the 8 targets covering the spectral range between [OII]A3727 
and [OIII]A5007. The lower panel shows the number of spectra used for each pixel, as the pixels interested by bright sky lines were not 
included. The positions of the undetected lines [NeIII]A3869 and HeI3889 are shown for reference. 



the 8 LSD sources whose spectra cover the [OII]A3727- 
[OIII]A5007 region. This spectrum was obtained as an av- 
erage of the original ones without any weighting, and 
reveals the average properties of the LSD sample as if 
they were observed simultaneously. Large [OIII] A5007/H/3 
and [OIII] A5007/ [Oil] A3727 ratios relative to local galax- 
ies are shown, and this is related to metallicity as ex- 
plained in section 14.61 The corresponding average mass is 
log(M./M0)=9.8O. 



6 ANALYSIS 

In the previous sections we have measured stellar mass, SFR, 
metallicity and size of our sample of LBGs. In this section 
we intend to use these parameters to derive information on 
the history and evolutionary state of these galaxies. We will 
show that mass-metallicity relation, gas masses and effective 
yields can provide important clues on the physical processes 
dominating the evolution of these objects and on the rela- 
tion between these starburst and similar objects at lower 
redshifts. 



6.1 The Mass-Metallicity relation 

Figure[5]shows the stellar mass-metallicity relation at 2~3.1, 
compared to the same relation as measured at l ower red- 
shifts. All the presented data have been scaled to a lChabrieil 
l|2003l ) IMF and use the same metallicity calibration. A 
strong, monotonic evolution of metallicity can be seen, i.e.. 



galaxies at 2~3.1 have metallicities ~6 times lower than 
galaxies of similar stellar mass in the local universe. It is 
worth noticing that this is not the evolution of individual 
galaxies, as discussed in iMaiolino et ahl l|2008l) . but this is 
the evolution of the average metallicity of the galaxies con- 
tributing to a significant fraction of the star-formation ac- 
tivity at their redshifts. The observed evolution implies that 
galaxies with relatively high stellar masses (log(M/M0)=9- 
11) and low metallicity are already in place at z >3, and 
this can be used to put strong constraints on the processes 
dominating galaxy formation. 

While stellar mass, based on integrated photometry, 
is representative of the full galaxy, metallicity is possibly 
dominated by the central, brightest regions. The presence of 
metallicity gradients could have some influence on the ob- 
served mass-metallicity relation. These aperture effects are 
present at any redshift: even at z~0, galaxy spectra from 
SDSS refer to the central few arcsecs of the galaxies. In 
most models, the central brightest part of the galaxies are 
also the most metal rich, therefore the use of total metal- 
licities for our LSD galaxies is expected to produce an even 
larger evolution. 



The effect of "downsizing" (jCowie et al.lll99"^ on chem- 
ical enrichment is expected to produce differential evolu- 
tion related to stellar mass. Stronger evolution for low-mass 
galaxies is observed from 2=0 to z=2.2 (see Figure [5]). The 
observed spread of the distribution and the uncertainties on 
the single points still make it impossible to see if such an 
effect is already in place between z=2.2 and z—S-A. Con- 
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Evol ution of the mass-metallicity relation from ^=0.07 llKewlev Ellisonll2008t) . to z=0.7 ||s avaglio et al. II2OO5I). 2= 2.2 

llErb et al.ll2006al) and 2=3-4 (AMAZE+LSD). All data have been calibrated to the same metallicity scale and IMF llChabrieJl2003) in 
order to make all the different results directly comparable. Turquoise empty dots show the AMAZE galaxies, blue solid dots the LSD 
galaxies. The solid square shows the "average" LSD galaxy, having average mass and composite spectrum (see Figure |4j. The lines show 
quadratic fits to the data, as described in the text. 



straints on this effect can be derived when the full AMAZE 

data sample will be presented. 

Using the same representation as in iMaiolino et al.l 

l|2008h . we fit the evolution of the mass-metallicity relation 
with a second order polynomial: 

Z = A [log(M,) - log(Mo)]' + Ko 

where A=-0.0864 and Mo and Kq are the free parameters of 
the fit. By using the LSD and AMAZE galaxies, we derive 
log(Mo)=12.28 and Kq = 8.69. The values of Mp and Kg for 
the sa mples at lower redshifts can be found in lMaiolino et all 
(|2008l ). and can be converted to the present system by sub- 
tracting log(1.7) to Mo. 

Several published models of galaxy formation (e.g., 
Ide Rossi et al.ll2007l : iKobavashi et al.1l2007f ) cannot account 
for such a strong evolution. The physical reason for this 
can be due to some inappropriate assumption, for example 
about feedback processes or merging h istory. When taken 
at face value, some o ther models (e.g., iBrooks et al ] |2007l : 
iTornatore et al.ll2007l l provide a better match with the ob- 
servations, but a meaningful comparison can only be ob- 



tained by taking into account all the selection effects and 
observational biases, and by comparing not only stellar mass 
and metallicity but also all the other relevant parameters, 
such as dynamical mass, angular momentum, gas fraction, 
SFR, morphology and size (see, for example, ICalura et al.l 

liboay 

In fact, it is important to emphasize that the galaxy 
samples used f o r Fig ure [5] change with redshift. In the 
iTremonti et al.l (|2004 ) work, the local SDSS galaxies un- 
der study constitute an almost complete census of the local 
star-forming galaxies, and the derived mass-metallicity re- 
lation shows the average properties of the sample. At high 
redshift, especially at z > 1, a number of effects must be 
considered. Only the most active galaxies are selected, and 
the observed metallicity refer to these objects. If, for ex- 
ample, less active, already formed and enriched galaxies are 
present at 2=3, they would not be present in our sample; the 
fraction of galaxies selected is likely to change with mass, 
and this is expected to introduce some systematic effect 
with stellar mass; the rest-frame UV selection misses extin- 
guished galaxies, therefore an increasing fraction of dusty. 
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Figure 6. SFR as a function of stellar mass for the galaxies in the LSD sample. This distribution is compared wi th two observed correla- 
tions at lower redshift (z~2), in magenta the UV-selected sample by Erb et al. (2006c), in green the SINS galaxies llForster Schreiber et al] 
|2009|) . The thick lines show the median of the correlations, while the color shaded regions are delimited by the 5%, 25%, 75% and 95% 
percentiles of each distribution. Compared to the sample at 2~2, LSD galaxies show a larger dispersion. The long-dashed lines show 
some values of star formation timescales, defined as t»=M«/SFR. LSD galaxies have SF timescales between 10^ and 10^ yr. 



metal rich galaxies are excluded from the mass-metallicity 
relation at high metallicities; the UGR color-color selection 
for the LBGs excludes galaxies with red (G-R) colors due 
to the presence of o lder stellar populations. For example, 
iHavashi et al.l l|2008l ) observed a sample of K-selected, star 
forming galaxies at 2~2. They found a mass-meta llicity re- 
lation which is similar to that in lErb et al.l (|2006al ) but with 
an offset of about 0.2 dex towards higher metallicities. The 
difference is likely to be due to target selection. Evolution in 
the mass-metallicity plane can also be produced by differ- 
ent SF histories of individual galaxies: if, for example, each 
galaxy experiences several bursts of star-formation, galaxies 
detected during the first burst will show low stellar masses 
and low metallicities, while galaxies with the same baryonic 
and total masses but detected during later bursts will show 
higher metallicities and stellar masses. This is not the case 
if galaxies have only one major burst, and all the galaxies 
are selected during the same phase. 

All these effects are present at any redshift but their 
importance changes significantly among the different sam- 
ples. For this reasons, the quantitative interpretation of the 



observed evolution requires great care. A meaningful com- 
parison between a sample of galaxies caught during a special 
phase of their life at «=3, with a more representative sam- 
ple of galaxies at z=Q can only be done when the selection 
effects are taken into account. This will be the subject of a 
future work. 



6.2 SFR and stellar mass 

Several studies of galaxies have found a correlation between 
stellar mass M* and SFR, i.e., more massive galaxies also 
have larger SFRs (see, for cxa,mplc, Scliiminovich ct al. 200?! 
for a UV-selected s ample. lElbaz et al.ll2007l for an optical se- 
lected sample, and lDaddi et al. I l2007l for a K-band selected 
sample). In contrast, local Ultra-luminous Infrared Galax- 
ies (ULIRG) do not follow the same relation but lie above 
it l|Elbaz et al.l |2007| ) . The slope of the correlation is not 
well constrained. While some authors (jPaddi et al.l |2007| ) 
find SFRs proportional to M., i.e., specific SFR (SSFR = 
SFR/M*) constant with mass, in most cases a much weaker 
dependence of SFR on M,is found (jNoeske et al.l l2007l : 
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I Zheng et all l2007l : iDrorv fc Alvare j I2OO8I '). Other authors 
find an evolution of the slope with redshift (|Dunne et alj 
l2008h . while in some samples the SFR seems to be more di- 
rectly relBted_tostellax surface density rather than stellar 
mass (|Franx et al.ll2008l ). 

The existence of such a relation and its slope are 
strongly affected by sample selection. For example, color- 
based selections tend to generate some proportionality be- 
tween M* and SFR because stronger SFRs are needed to 
significantly affect the colors of more massive galaxies. For 
UV-selected galaxies with SFRs estimat ed from th e optica l 
emission lines, similar to LSD, at 2:~2 lErb et al.1 (|2006al l 
have found a positive, albeit weak, correlation between M* 
and SFR by removing the galaxies which have dynamical 
mass much hi gher than the stellar rnass. I n contrast, the 
SINS sample jForster Schreiber et af] 12009 1') shows a much 
stronger dependence of SFR on M, with slope consistent 
with one, in agreement w ith the corre l ation found for the 
K-band selected sample in iDaddi et al.1 (|2007l ). 

To produce such a correlation, galaxies must spend a 
large portion of their life forming stars at a characteristic 
fractional level, instead of forming most of their stellar mass 
in a few large bursts. As a consequence, this correlation 
is often interpreted as an evidence that most of the star- 
formation activity is not related to short, intense burst but 
rather is associated to a more continuos activity at a lower 
level. This i s confirmed, at 2:^ 2.2, by the fact that the cur- 
rent SFR in lErb et all (|2006d ) is similar to the past average. 
This ongoing activity could be rel ated to cold accretion of 
gas followed by d isk instabilities ijPekel fc Birnboiml l2006l : 
iDekel et al.ll2009l ). The regular morphologies of the galaxies 
domin ating the star formation activity at intermediate red- 
shifts ()BeU et al.ll2005l : IWoff et al.ll2005ll and the (generally) 
regul a r dynamics observed at z~2 jForste r Schreiber et al.l 



20061 : iGenzel et all l2008l : IShapiro et al.1 12008: Stark ct al 



20081 ) support this interpretation. 



It is interesting to study if a correlation is present in our 
sample, and Figure [6] shows the relation between M, and 
SFR in the LSD sample. Large uncertainties are present, in 
particular on the SFR that is very sensitive to the poorly 
constrained amount of dust extinction. The LSD galaxies 
appear to have SFRs and M* similar to the lower-redshift 
galaxies in the samples mentioned above, but the small num- 
ber of data points and the large intrinsic errors makes it im- 
possible to obtain robust conclusions about the existence of 
a correlation. 

It is important to study if selection effects could hide a 
strong correlation between M* and SFR. The rest-frame UV 
selection explains why no galaxy below a given SFR thresh- 
old is selected. On the contrary, the lack of active, massive 
could be due to several effects: 1- high-mass galaxies selected 
as LBGs with high SFRs could be rare objects, and therefore 
none of them is included in this sample of 10 objects; 2- it is 
possible that the dust content of the galaxies increases with 
age and metallicity, i.e., with stellar mass. In this case mas- 
sive galaxies with large SFRs could be too red or too faint 
in the UV to be selected as LBGs, and could be present in 
other catalogs as sub-mm galaxies or ULIRGs. Fo r example, 
this is what is observed bv lNesvadba et aP l|2007h in the the 
archetypal sub-mm galaxy SBS J14011-I-0252, which is both 
metal rich, actively star-forming, and dusty; 3- the SF ac- 
tivity of each galaxy could be a decreasing function of time. 
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Figure 7. Mass of gas as a function of the stellar mass for 
the LSD galaxies. The dotted lines shows equal masses. All the 
galaxies have similar amount of gas, which is not related to stellar 
mass. 



with a strong burst followed by a rapid decrease. In this case 
massive galaxies would also have low SFRs. 

Albeit these limitations, our data do not support the 
existence of an almost direct proportion ality (SFR ~M«^'° ) 
bet ween SFR and M« similar to t hat in lDaddi et al.1 (|2007l ) 
and lForster Schreiber et al.l (|2009l ) . In contrast, we find bet- 
ter agreement with a fiatter relation (SFR ~M«° '') as found 
bv lErb et all f(2006a ). This implies the presence of larger 
SSFRs in lower meiss galaxies. As shown by the dashed 
lines in Fig. (S] the same effect can be expressed by the 
star-formation timescale, defined as the time to form 
the present amount of stars M« at the present level of SF, 
t*=M*/SFR=l/SSFR. In the LSD sample, t* changes sys- 
tematically with Af«: low mass galaxies have ^ lO^yr, 
while higher mass galaxies have ~ lO^yr. This is an in- 
dication that, in relative terms, the current episode of star 
formation is more important in lower-mass systems. While 
high-mass galaxies have SFRs and SSFRs similar to what 
observed at lower redshift, low-mass galaxies tend to be 
relatively more active. This interpretation implies that less 
massive galaxies are also younger. An accurate test of this 
prediction is not easy, as it is difficult to measure the age of 
high-redshift galaxies. Nevertheless, it is consistent with the 
broad correlation between stellar mass and age described in 
sect. 1331 



6.3 The mass in gas 

In local star-forming galaxies, the surface density of star 
formation F^sfr is related to the surface density of gas Togas'- 

This is t he Schmidt-Kennicutt law (see, for example, 
lKennicut t| |1998), which seem s to be already in place at high 
redshifs (|Bouche et al.ll2007l ). It is interesting to note that 
downsizing is a natural consequence of the Schm i dt law 
and of the other scaling relations, as shown bv lErbl (|2008l ). 
If infalls and outfiows can be neglected on the timescale of 
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a starburst, or if they are proportional to the SFR itself, 
the Schmidt law produces a decreasing SFR whose initial, 
maximum value scales with the mass M and the size r 



as A/i 



and the typical timescales as M 



If we assume that galaxies follow the same mass-size 
rel ation observed in the local massive, early-type galaxies 
bv lShen et all (|2003h . we obtain that the typical timescale 
decreases with mass as M""'^ and, correspondingly, the 
magnitude of the initial starburst increases as M^-^ The 
typical timescales computed by lErbl l|2008l ') are ~0.3 Gyrs 
for M = lO^^Mo, and ~4 Gyr for M = IO^Mq. 

The Schmidt-Kennicutt law is often used to estimate 
gas density and gas mass sta rting from the observed SFR 
and galaxy size (for example, lErb et al.ir2006al : lErbll2008h . 
The slope n of the relation is still unc e rtain. The "classic" 
value of n is 1.4, given by iKennicutO (Il998l ). Star forma- 
tion due to lar ge-scale instabilit ies of the disk naturally pro- 
duces n= 1.5 (|ElmegreeDll2002l ). an d such a value is sup - 



ported by more recent observations (iKennicutt et al.1 



2007 ) 



on small scales. In high redshift studies. iBouche et al.l (|2007l ) 
found a slightly larger value, n=1.7, a value which depends 
on the adopted value of the CO/H2 ratio. Such a value 
could be more appropriate for our LBGs, nevertheless for 
consistency with previous works we adopt the relation in 
iKenniciitt (|l998l l. deriving the following equations for the 
gas surf ace d e nsity S t, g a and the gas mass Mgas, valid when 
the lChabrieJ (|2003l ) IMF is used: 



Egas{MQ/pc) = 254 



SFR 



Mo/yr 



kpc 



SFR 



(1) 



(2) 



Using the values in iBouche et all (|2007l ) would result in 
gas fractions about 40% lower, which however do not affect 
the conclusions of this work. 

The largest uncertainties in these equations are related 
to the half-light radius r, which has large errors, as explained 
above. This is important for Egas, which has a strong depen- 
dence on galaxy size. In contrast, the total gas mass Mgas 
is much less sensitive to galaxy size as it depends on r"'^® 
and, as a consequence, the uncertainties on this quantity are 
smaller. 

The observed SFR densities are between 0.5 and 
6Mp) /yr/ kpc'^, similar to what observed both in the loca l 
universe (|Kennicuttlll998l ') and at z=2.2 (lErb et al.ll2006d ). 
The corresponding gas densities range between 300 and 2000 
M (7)/pc^ , similar to the valu es derived for LBGs a nd ULIRGs 
bv lCoppin et al.l (|2007f ) and lTacconi et all (|2006l '). As we are 
sampling the central, most active part of the galaxies, these 
numbers must be considered as the maximum surface den- 
sities of gas. 

Figure [7] shows the resulting Mgaa as a function stellar 
mass. It is evident that gas mass does not correlate with 
M« and that Mgas covers a much smaller range than Af, . 
While the largest galaxy has a stellar mass 24 times larger 
than the smallest, the most gas-rich galaxy has only 5 
times more gas than the poorest. This result could be 
partly due to a combination of selection effects and spatial 
resolution. As all the galaxies have similar SFRs, the range 
in gas mass can be dominated by the range in intrinsic 



dimensions. Our resolution, 0.2" FWHM, corresponding to 
about 1.4 kpc, is similar to the intrinsic galaxy size. As 
a consequence, it is possible that we are overestimating 
the size, and therefore the gas mass, of the smallest galaxies. 

Interesting information on the physical properties of the 
LSD galaxies can be obtained by comparing the different 
timescales. We have already introduced the stellar timescale 
t,, which measures the time to create the observed amount of 
stars at the current SFR. The typical dynamical timescale is 
tdyn = 2nR'^^^ /(GMy^^ , where R is the typical radius and 
M the typical total ma ss. Any starburst activity is expected 
to last at least t^yn l|Dekel fc BirnboimI [2009 1. The LSD 
galaxies have tdyn between 0.4 and 1.0x10^ yr. This nar- 
row range of values is not surprising, because if size depends 
on mass as in local early-type galaxies, dynamical times are 
almost independent of mass. A third important timescale is 
the gas exhaustion time tex, i.e., the time to exhaust all the 
available gas at the current rate of SF, tex ~ Mgas/ SFR. 
For the LSD galaxies, this spans a narrow range, between 3 
and 6x10* yr. 

We have already seen that Mgas/M,, and therefore 
t*/tex, varies systematically with M*, meaning that the 
smallest galaxies are more gas rich and are forming a larger 
fraction of stars. The ratio between tdyn and tex ranges 
between 6 and 4, with no clear dependence on M*. This 
implies that the galaxies are so gas rich that the current 
level of star formation can be sustained for several tdyn- We 
cannot be sure that the SF activity will proceed for so long, 
but this is an indication that we are dealing with major 
episodes of star formation in galaxies which have accreted 
large amounts of gas, of the order of 10^" M© (see Fig.[7|. In 
these galaxies star formation follows mass assembly. 



6.4 The fraction of gas 

The gas mass Mgas can be compared with the total baryonic 
mass Mh to estimate the fraction of mass in gas, an indica- 
tion of the evolutionary stage of the galaxy. The results are 
shown in Figure [S] where the gas fraction is plotted vs. the 
stellar mass. It should be noted that stellar and gas masses 
are not measured in the same apertures. When HST images 
for these objects are available (see sect. l3.2| l. we will be able 
to quantify this effect and apply a correction. Despite the 
large uncertainties, a clear correlation is seen between stellar 
mass and gas fraction. This is similar to what is observed at 
low-redshifts, where a tight cor relation exists b etween gas 
fracti on and rotation velocitv (Dalcantonl |2007| ) or stellar 
mass (|Bellll2003l: iKannappanlbooi T The same effect is also 
observed bv lErb et al.1 (|2006al ') and iReddv et al] (|2006l ) in 
their samples at z~2. The relation between stellar mass and 
gas fraction in local galaxies is o ften used to constrain the 
models of galaxy formation (e.g., [Somerville et al.ll2008l ). 

The observed correlation is at least partly due to the se- 
lection effects of our target sample, shown as a grey-shaded 
area. As explained above, our sample is selected according 
to the SFR, while stellar mass can assume a wide range of 
values. The existence of this selection effect makes it impos- 
sible to test at z= 3 the results of seve ral models of galaxy 
formation, such as lBrooks et al.l l|2007h and lMouchine et al.l 
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Figure 8. Fraction of baryonic mass in gas as a function of stellar 
mass. The blue solid dots are the LSD galaxies. The grey-shaded 
area represents galaxies that cannot be included in our sample 
because the SFR is lower than our limit (~15 M-q/jt). The dashed 
lines show the expected dependence of gas fraction on stellar mass 
for galaxies having SFR constant (green) or proportional to the 
stellar mass (red). 



l)2008t ). These models predict that a strong correlation ex- 
ists between fgae and M« at any redshift, with more massive 
galaxies having a lower conte nt in gas, very simil ar to what 
is actually observed (see also ICalura et al]l2008l ). In other 
words, the shaded region in Fig. [8] could be empty for phys- 
ical reasons. 

The expected gas fraction as a function of stellar mass 
can be derived once the dependence of dimensions and 
SFR on M-t is known. We compute this gas density using 
the mass-s i ze rel ation for massive, late- type galaxies in 
IShen et all (|2003l ). r oc M° *, and assuming that the SFR 
of the galaxies does not depend on stellar mass. As size 
are expected to change slowly with mass, these results 
are largely independent of the evolution of the mass-size 
relation, which is largely unknown. The result is plotted 
as a short-dashed, green line in Figure [51 In contrast, 
assuming that the SFR is proportional to stellar mass 
as observed in the local universe (see section 14. 5p . we 
derive the long-dashed, red line in the same figure. The 
normalization of these two lines can vary as a function of 
the adopted SFR and galax;y size, but the shape of the 
curves show that the gas fraction measured in our target 
galaxies is that expected for a SFR-limited sample. 



6.5 The effective yields 

Important hints on the physical processes shaping the mass- 
metallicity relation can be obtained by considering the effec- 
tive yields, i.e., the amount of metals produc ed and retained 
in the ISM per unit mass of formed stars l|Garnett|[2003 : 



Figure 9. Effective yields of the LSD sample (solid blue dots), as 
a function of stellar mass, compared with the solar value (dotted 
line) . Light and dark error bars show the contribution to the errors 
on yields due to the uncertainties on gas fractio n and metallicity, 
respectively. The magenta stars are the results bv lErb et al .l2006al 
at 2=2.2. The black solid line are the local results by Tremonti 
et al. (2004). The green dashed line corresponds, as in Figure [S] 
to galaxies with constant SFR. The blue, long-dashed line is the 
model with infalls and outflows proportional to the SFR described 
in the text. The labels near the top of the figure show the values 
of the average metallicity for galaxies of a given stellar mass as 
derived from the mass-metallicity relation. 



lDalcantonll2007h . In a closed-box model, metallicity Z (i.e., 
the mass in metals divided by the mass of the interstellar 
gas) can be written as a simple function of the gas fraction 

fgas '■ 



(3) 



where y is the true stellar yield, i.e., the ratio between the 
amount of metals produced and returned to the ISM and 
the mass of stars. This is the solution of a differential equa- 
tion (see .Edmunda i,1990) valid in a closed-box case, with in- 
stantaneous recycling, instantaneous mixing, and low metal- 
licities. The solar yield j/0, i.e., the fractional contribution 
of metals to the solar mass, is often used as reference and 
is 0.0126 (jAsplund e t al. 2004*1. By inverting this equation, 
and using observed quantities for Z and fgas, we obtain the 
effective yields: 

= Z/ln{l/fgas) 



(4) 



The measured values of j/e/ / could differ from the true 
stellar yields y if some of the assumptions used to derive 
eq. |3]do not hold, in particular if the system is not a closed 
box. Outflows of very enriched material, such as the ejecta 
from SNe, extract metals, decreasing Z, while infalls of 
pristine gas can increase fgas and decrease Z. It is well 
known that both effects must be present at some level. Gas 
infall, either in merging episodes or in cold gas accretion, 
is needed to supply galaxies with the amount of matter 
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we observed in stars today, while starbursts are known to 

J iroduce SN driven winds of the order of 0.1 MQyr~^ kpc~^ 
Lehnert fc Heckmanlll996al ). The presence of both effects 
does not imply that both contribute significantly to the 
chemical evolution of the systems and, in particular, to the 
observed yields. In all case s, yeff must be sm alle r than y, as 
discussed, for example, bv lEdmun ii l|l990l ) and iDalcantonI 
l|2007l ). This is why it is interesting to measure the effective 
yields, because it can measure how a system is far from 
being a closed-box. 



O.B 



In the local universe, iTremontT et all (|2004l ') have found 
a significant dependence of yeff on mass, i.e., lower mass 
galaxies have lower j/e//. The reduction of y^ff below the 
solar value in galaxies with low stellar mass is usually 
described as a consequence of outfiows. If outfiows are the 
main effect shaping the mass-metallicity relation, then the 
effective yields are expected to increase with stellar mass be- 
cause lower mass galaxies, having a shallower potential well, 
have lost a larger fraction of metals into the intergalactic 
medium. Above a certain stellar mass, the galaxy potential 
well is too deep for the SNe to eject significant fraction of 
enriched gas, and J/e//~ y is expected in massive galaxies. 

Figure |§] shows the values of yef / for the LSD galax- 
ies. The errors on each single point are large, as they reflect 
the large uncertainties on both gas fractions and metallici- 
ties. Systematic effects on dimensions could also be present, 
as discussed above. Nevertheless, the situation seems pretty 
clear: rather than increasing with stellar mass toward solar 
values, in the LSD sample y^fi is found to decrease with 
stellar mass, starting from vp, at the low mass end . The dif- 
ferences with the results bv lTremonti et al] (|2004h are strik- 
ing. This is not surprising because it means that the physical 
processes dominating the SDSS galaxies in the local universe 
are not the same as in LBGs at z >3. In contrast, our resu lt 
is very similar to what has been found bv lErb et all (|2006al ). 
Our and Erb's samples share the same dependence of yefj 
with stellar mass, although LSD galaxies have lower y^ff 
because of thei r lowe r metallicity. Our results ar e also simi- 
lar to those of iWeiner et al. (200^) at 2=1.4 and lLaw et al.l 
at 2=2.5, who detected the presence of stronger galac- 
tic winds in more massive galaxies. 

To interpret this result we have to take into account how 
the samples were chosen and the scaling relations. From the 
expected gas fraction in Figure |5] and the mass-metallicity 
relation in Figure O we can predict the j/e// of the LSD 
sample. This is a real prediction, i.e., once you define SFRs, 
metallicity and dimensions, there are no free parameters to 
compute yields. We plot this prediction in Figure |9l The 
observed behavior is in agreement with what is expected 
for a SFR-limited sample. In other words, the dependence 
of the yields on stellar mass is a direct consequence of the 
mass-metallicity relation and of the other scaling relations. 

The presence of yields decreasing with mass implies that 
the classical outflows, whose specific power decreases with 
increasing galaxy mass, do not apply to the LBG at 3.1. 
The mass-metallicity relation at z >3 must have a different 
origin. 

The dependence of j/e// on stellar mass is a secondary 
effect due to the relation between stellar mass and gas frac- 
tion. For this reason, the physical meaning of these findings 
is better understood by plotting y^f / directly as a function of 
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Figure 10. Metallicity [top] and effective yields (bottom) as a 
function of gas fraction. This is plotted decreasing rightwards, to 
have galaxies evolving from left to right. The blue solid dots are 
the LSD galaxy sample, the magenta stars the lErb et al. I ll2006al ') 
data. The black dashed line shows the expectations from a closed- 
b ox model. T he solid lines show the results of applying the model 
in lErbl | |2008|) with values of fi and /o in the labels. The red line 
shows a model with pure outflows, with ejected mass 6 times the 
SFR. The blue line is a model with both infall (/i=2.5) and out- 
flow (/o=l). The blue dotted line is this same model but without 
outflow {/i=2.5 and /o=0). 



gas fraction. This is done in Figure [TOl The relation between 
these two quantities and the effect of infall s and outflows are 
described by several authors. In particular. IDalcantonI (|2007| ) 
shows that outflows of very enriched materials can be very 
effective in reducing yeff, especially for gas-rich systems. 
Gas-poor system can change their y^ff by a factor of two, 
as observed in the LSD galaxies, if they accrete a significant 
amount of gas, larger than the internal gas. 

For gas rich systems, j/e// do es not depend on fgas and 
it is equal to the stellar yields y (|Dalcantonll2007l ), thus al- 
lowing us to measure this quantity, which in principle could 
be different from y©. Metallicity can infiuence nucleosynthe- 
sis in different ways. The most important effect is expected 
to be mass-loss, which has a strong dependence on metal- 
licity. Massive stars with stronger stellar winds lose a larger 
fraction of He and C that would otherwise be converted to O. 
As a consequence, the oxygen yield we measure is expected 
to decrease with metallicity. The expected effect is not large, 
up to a factor of 2, but in principle measurable with our data 
(jWooslev fc WeaveJll995f : iMarigollioOll : lGarnettll20o3 ). AU 
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the LSD galaxies with fgas> 0.7 have y^ff fully consistent 
with j/0, showing that any metallicity dependence of y is not 
large even eX Z ^ O.IZ©. 



iMatteuccil (j200ll ) (see also iMatteuccil 120081 ;) and lErbl 
('2008') introduce a model in which galaxies have infalls and 
outflows proportional to the SFR. The instant recycling and 
mixing approximations are used, infalling gas is assumed to 
have no metallicity, while the metallicity of the outflowing 
gas is considered to be the same as the ISM. This model can 
be used to reproduce j/e// by changing two free parameters, 
i.e., the amount of infall ji and outflow /□ in unit of the SFR 
of the galaxy. Two more parameters can be varied, if nec- 
essary, i.e, the true yield y, and the fraction a of mass still 
locked in stars. Implicitly, the physical properties of these 
outflows differ significantly with what is invoked to explain 
the yields in t he local universe: given the correlation betw een 
SFR and M. (|Schiminovich et al.|[2007l : lErb et al.ll2006d '). at 
2~2 these outflows are exp ected to increase , rat her than de - 
crease, with stellar mass. ,Erb et al.1 (1200681 ) and lErbl (120081 ) 
find that their data at z=1.2 can be adequately reproduced 
by a super-solar true yield (t/=1.5y0) and significant infalls 
and outfiows, with the best-fitting parameters of /i=2.2 and 
/<,=1.3. 

We have applied this model to LSD, reproducing the 
dependence of ?/e// on jgas varying /i and fo- Infalls are 
effective in changing y^// for the gas-poor galaxies, leaving 
gas-rich galaxies much more unaffected. Outfiows are more 
efficient in reducing ye// in gas-rich galaxies, a nd m u ch les s 
effective for gas-poor gala:xies. In contrast to lErbl (|2008t ). 
we have fixed the values of y to y^, as explained above. 
The parameter a is a slowly decreasing function of time, as 
more stars leave the main sequence, but can vary in quite a 
small range. It is expected to be ~0.93 at t=10^yr, ~0.84 
at t=10*yr, ~0.76 at t=10'^yr, and approachin g 0.60-0.65 
after an Hubble time (|Bruzual fc Charlot|[2Cio3 ). Given the 
ages of our galaxies, we have used a=0.8, and the results are 
not very sensitive to the value chosen. The results are shown 
in Figure [TUl Outfiows without infalls can explain the ob- 
served 2/e/ / , but in this case large values of /□ between 4 and 
8 are needed. The Figure shows the best-fitting pure out- 
fiow model, having /o=6. In other words, if pure outfiows 
are responsible for the reduction of j/e//, the most active 
galaxies must eject into the IGM masses of the order of 400- 
8OOM0 /yr, while converting IOOMq /yr into stars. This must 
happen in galaxies as massive as 3 • lO^^Mo, i.e., galaxies 
that , in the local univers e, show solar yields. As a compari- 
son, IWeiner et al.l (|2008| ) found outflowing masses of about 
20 Mo/yr in their sample of starburst galaxies at 2=1.4, 
whose SFRs are 10-100 M© /yr, similar to the present sam- 
ple. Pure infalls can explain the yields with lower amount 
of flowing gas. In this case, the best fitting value is /i=2.5, 
with an acceptable range between /i=1.5 and /i=3.5. Once 
an infall with /i ~ 2 is present, it dominates the behavior 
of ye// and jo can assume any value between and 7. The 
best-fitting model with fi— 2.5 and fo~l is also plotted in 
Fig. |5]as a function of stellar mass, providing a good fit. 



7 SUMMARY AND CONCLUSIONS: A 
PICTURE OF LEGS 

We have obtained deep, spatially resolved, near-IR spec- 
troscopy of a complete sample of LBGs at z~3, selected 
only to be near a bright foreground star needed to drive the 
adaptive optics system. These gal axies are expected t o give 
an unbiased representation of the lSteidel et al.l (|2003l ) sam- 
ple of LBGs. The detected optical lines give information on 
many aspects of these galaxies, and in this paper we have 
analyzed the properties related to metallicity, stellar and gas 
mass, and gas fraction. 

A strong evolution of the mass-metallicity relation is 
found, i.e., LBGs at 2:~3 have metallicities between 0.1 and 
0.5 solar. Assuming that the Schmidt-Kennicutt law is al- 
ready in place at 3, we derive the properties of the gas. 
Large masses of gas (~ IO^'^Mq) are present in these galax- 
ies, as a result of recent episodes of gas accretion. The gas 
fraction is found to be anti-correlated with stellar mass, and 
this is possibly due to the UV selection of the sample. These 
results on gas mass and metallicity imply that, at these red- 
shifts, LBGs are already massive but still metal-poor. This 
means that mass assembling is occouring before star forma- 
tion, as observed in carly-typc galaxies at lower redshifts 
fe.g.. lSaracco et al.ll200ll2005l : lOaddi et al.ll2005D Once the 
selection effects are accurately considered, this can be used 
to put constraints on the models of galaxy formation and 
evolution. 

The time scale to produce the observed stars at the 
current level of star formation varies from ~ lO^yr for the 
smallest systems to ~ lO^yr for the galaxies with the largest 
stellar mass. In contrast, the timescale to exhaust the avail- 
able gas is less variable, and is comparable with the dynam- 
ical times of these systems. 

The effective yields decrease with stellar mass and have 
solar values at the low-mass end. This is in contrast with 
what is observed in the local universe, and can be ascribed 
to recent accretion of metal poor gas at a rate of the order 
of the current SFR. Stellar yields have solar values even at 
Z^O.IZq. 

All these evidences show that the LBGs of our sam- 
ple can be described as starbursts on several dynamical 
timescales following major events of gas accretion, due 
either to gas infall or merging with gas-rich, metal poor 
galaxies. Low mass galaxies can be described as expe- 
riencing their first major burst of star formation. More 
massive galaxies are more gas poor and have lower yields, 
as is expected when geis infall ignites a new episode of 
star formation in an older galaxy. Outflows are probably 
present but do not dominate the chemical evolution of these 
systems, that is more easily driven by infalls. 

These galaxies will soon be observed by HST at optical 
and near-IR bands. These data will complement our near-IR, 
spatially resolved spectroscopy in oder to study the morpho- 
logical and dynamical properties of these galaxies and obtain 
a complete picture of these LBGs. 

Acknowledgments We thanks F. Eisenhauer and A. 
Modigliani for support during data reduction, and B. Ken- 
nicutt and T. Nagao for very useful discussions. We also 
thanks the staffs of ESO, Spitzer and TNG for excellent 
service observing, and E. Daddi for having provided data in 



18 F. Mannucci et al. 



tabular form. This work was partially supported by the Ital- 
ian Space Agency through contract ASI-INAF 1/016/07/0, 
by INAF CRAM 1.06.09.10, and by NASA (Spitzer) grant 
number 1343503. 



REFERENCES 

Argcnce, B. & Larnarcillc, F. 2008, ArXiv e-prints 
Asplund, M., Grevesse, N., Sauval, A. J., Allende Prieto, 

C, & Kiselman, D. 2004, A&A, 417, 751 
Bafifa, C, Comoretto, G., Gennari, S., et al. 2001, A&A, 

378, 722 
Bell, E. F. 2003, ApJ, 586, 794 

Bell, E. F., Papovich, C, Wolf, C, et al. 2005, ApJ, 625, 
23 

Bolzonella, M., Miralles, J.-M., & Pello, R. 2000, A&A, 
363, 476 

Bonnet, H., Abuter, R., Baker, A., et al. 2004, The Mes- 
senger, 117, 17 

Bouche, N., Cresci, G., Davies, R., et al. 2007, ApJ, 671, 
303 

Brinchmann, J., Pettini, M., & Chariot, S. 2008, MNRAS, 
385, 769 

Brooks, A. M., Governato, F., Booth, C. M., et al. 2007, 
ApJL, 655, L17 

Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 

Calura, F., Jimenez, R., Pantcr, B., Matteucci, F., & Heav- 
ens, A. F. 2008, ApJ, 682, 252 

Calura, F., Pipino, A., Chiappini, C, et al. 2009, ArXiv 
e-prints 

Calzetti, D., Armus, L., Bohhn, R. C, et al. 2000, ApJ, 
533, 682 

CardelU, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 
345, 245 

Chabrier, G. 2003, PASP, 115, 763 

Chapman, S. C, Blain, A. W., Small, I., & Ivison, R. J. 
2005, ApJ, 622, 772 

Conselice, C. J., Bundy, K., Trujillo, I., et al. 2007, MN- 
RAS, 381, 962 

Coppin, K. E. K., Swinbank, A. M., Neri, R., et al. 2007, 

ApJ, 665, 936 
Cowie, L. L. & Barger, A. J. 2008, ApJ, 686, 72 
Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, 

AJ, 112, 839 

Crosci, G., Davies, R. I., Baker, A. J., & Lehnert, M. D. 

2005, A&A, 438, 757 
Cresci, G., Hicks, E. K. S., Genzel, R., et al. 2009, ArXiv 

e-prints 

Daddi, E., Dickinson, M., Morrison, G., ct al. 2007, ApJ, 
670, 156 

Daddi, E., Renzini, A., Pirzkal, N., et al. 2005, ApJ, 626, 
680 

Dalcanton, J. J. 2007, ApJ, 658, 941 

Davies, R. I. 2007, MNRAS, 375, 1099 

de Rossi, M. E., Tissera, P. B., & Scannapieco, C. 2007, 

MNRAS, 374, 323 
Dekel, A. & Birnboim, Y. 2006, MNRAS, 368, 2 
Dekel, A., Sari, R., & Ccvorino, D. 2009, ArXiv e-prints 
Drory, N. & Alvarez, M. 2008, ApJ, 680, 41 
Dunne, L., Ivison, R. J., Maddox, S., et al. 2008, ArXiv 

e-prints 



Edmunds, M. G. 1990, MNRAS, 246, 678 

Eisenhauer, F., Abuter, R., Bickert, K., et al. 2003, in So- 
ciety of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, Vol. 4841, Society of Photo-Optical In- 
strumentation Engineers (SPIE) Conference Series, ed. 
M. lye & A. F. M. Moorwood, 1548-1561 

Elbaz, D., Daddi, E., Le Borgne, D., et al. 2007, A&A, 468, 
33 

Ellison, S. L., Patton, D. R., Simard, L., & McConnachie, 

A. W. 2008, ApJL, 672, L107 
Elmegreen, B. G. 2002, ApJ, 577, 206 
Erb, D. K. 2008, ApJ, 674, 151 

Erb, D. K., Shapley, A. E., Pettini, M., et al. 2006a, ApJ, 
644, 813 

Erb, D. K., Steidel, C. C, Shapley, A. E., et al. 2006b, 

ApJ, 647, 128 

Erb, D. K., Steidel, C. C, Shapley, A. E., et al. 2006c, ApJ, 

646, 107 

Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 
154, 10 

Forstcr Schreiber, N. M., Genzel, R., Bouche, N., et al. 

2009, ArXiv e-prints 
Forster Schreiber, N. M., Genzel, R., Lehnert, M. D., et al. 

2006, ApJ, 645, 1062 
Franx, M., Labbe, I., Rudnick, G., et al. 2003, ApJL, 587, 

L79 

Franx, M., van Dokkum, P. G., Schreiber, N. M. F., et al. 

2008, ApJ, 688, 770 
Frye, B., Broadhurst, T., & Bem'tez, N. 2002, ApJ, 568, 

558 

Gallazzi, A., Chariot, S., Brinchmann, J., & White, 

S. D. M. 2006, MNRAS, 370, 1106 
Garnett, D. R. 2002, ApJ, 581, 1019 
Gavazzi, G. & Scodeggio, M. 1996, A&A, 312, L29 
Genzel, R., Burkert, A., Bouche, N., et al. 2008, ApJ, 687, 

59 

Halliday, C, Daddi, E., Cimatti, A., et al. 2008, A&A, 479, 
417 

Hayashi, M., Motohara, K., Shimasaku, K., et al. 2008, 
ArXiv e-prints 

Hopkins, A. M. & Beacom, J. F. 2006, ApJ, 651, 142 

Kannappan, S. J. 2004, ApJL, 611, L89 

Kennicutt, Jr., R. C. 1998, ARAA, 36, 189 

Kennicutt, Jr., R. C, Calzetti, D., Walter, F., et al. 2007, 
ApJ, 671, 333 

Kcwlcy, L. J. & Ellison, S. L. 2008, ApJ, 681, 1183 

Kobayashi, C, Springel, V., & White, S. D. M. 2007, MN- 
RAS, 376, 1465 

Kobulnicky, H. A. & Koo, D. C. 2000, ApJ, 545, 712 

Koopmans, L. V. E. & Treu, T. 2003, ArXiv Astrophysics 
e-prints 

Koppen, J., Weidner, C, & Kroupa, P. 2007, MNRAS, 375, 

673 

Lamareille, F., Brinchmann, J., Contini, T., ct al. 2008, 
ArXiv e-prints 

Law, D. R., Steidel, C. C, Erb, D. K., et al. 2007, ApJ, 
656, 1 

Lee, H., Skillman, E. D., Cannon, J. M., et al. 2006, ApJ, 

647, 970 

Lehnert, M. D. & Heckman, T. M. 1996a, ApJ, 462, 651 
Lehnert, M. D. & Heckman, T. M. 1996b, ApJ, 472, 546 



LSD: Mass, metalUcity and gas at z^S.l 19 



Lequeux, J., Peimbcrt, M., Rayo, J. F., Serrano, A., & 

Torres- Peimbert, S. 1979, A&A, 80, 155 
Liu, X., Shapley, A. E., Coil, A. L., Brinchmann, J., & Ma, 

C.-P. 2008, ApJ, 678, 758 
Maier, C, Lilly, S. J., CaroUo, C. M., at al. 2006, ApJ, 639, 

858 

Maiolino, R., Nagao, T., Grazian, A., at al. 2008, A&A, 
488, 463 

Mannucci, F., Buttery, H., Maiolino, R., Marconi, A., & 

Pozzetti, L. 2007, A&A, 461, 423 
Mannucci, P., Pozzetti, L., Thompson, D., at al. 2002, MN- 

RAS, 329, L57 
Maraston, C. 2005, MNRAS, 362, 799 
Marigo, P. 2001, A&A, 370, 194 

Matteucci, F. 2001, Astrophysics and Space Science Li- 
brary, Vol. 253, The chemical evolution of the Galaxy 
(Kluwer) 

Matteucci, F. 2008, ArXiv c-prints 

Michcl-Dansac, L., Lambas, D. G., Alonso, M. S., & Tis- 

sera, P. 2008, MNRAS, 386, L82 
Modighani, A., Hummel, W., Abuter, R., et al. 2007, ArXiv 

Astrophysics e-prints 
Mori, M., Umamura, M., & Farrara, A. 2004, ApJL, 613, 

L97 

Mouchine, M., Gibson, B. K., Randa, A., & Kawata, D. 

2008, ArXiv e-prints 
Moustakas, J., Koruiicutt, Jr., R. C., & Tremonti, C. A. 

2006, ApJ, 642, 775 
Nagao, T., Maiolino, R., & Marconi, A. 2006, A&A, 459, 

85 

Ncsvadba, N. P. H., Lehnart, M. D., Eisanhauar, F., et al. 

2006, ApJ, 650, 661 
Ncsvadba, N. P. H., Lehnart, M. D., Ganzal, R., at al. 2007, 

ApJ, 657, 725 

Noeske, K. G., Faber, S. M., Weiner, B. J., at al. 2007, 

ApJL, 660, L47 
Pantcr, B., Jimenez, R., Heavens, A. F., & Chariot, S. 2008, 

MNRAS, 391, 1117 
Perez-Gonzalez, P. G., Gil de Paz, A., Zamorano, J., et al. 

2003, MNRAS, 338, 525 
Perez-Montero, E., Contini, T., Lamarailla, F., at al. 2008, 

ArXiv e-prints 

Pattini, M., Rix, S. A., Staidal, C. C., et al. 2002, ApJ, 569, 
742 

Pettini, M., Shapley, A. E., Stcidel, C. C, et al. 2001, ApJ, 
554, 981 

Pilyugin, L. S., Vi'lchez, J. M., & Contini, T. 2004, A&A, 
425, 849 

Pozzetti, L., Bolzonalla, M., Lamarailla, F., at al. 2007, 

A&A, 474, 443 
Reddy, N. A., Staidal, C. C., Fadda, D., at al. 2006, ApJ, 

644, 792 

Rodriguas, M., Hammer, F., Flores, H., et al. 2008, A&A, 
492, 371 

Saracco, P., Longhetti, M., Severgnini, P., et al. 2005, MN- 
RAS, 357, L40 

Saracco, P., Longhetti, M., Severgnini, P., et al. 2003, A&A, 
398, 127 

Savaglio, S., Glazebrook, K., La Borgna, D., at al. 2005, 
ApJ, 635, 260 

Scannapiaco, C., Tissara, P. B., White, S. D. M., & 
Springal, V. 2008, MNRAS, 389, 1137 



Schiminovich, D., Wydcr, T. K., Martin, D. C, ct al. 2007, 

ApJS, 173, 315 
Shapiro, K. L., Genzel, R., Forster Schreiber, N. M., et al. 

2008, ApJ, 682, 231 
Shapley, A. E., Coil, A. L., Ma, C.-R, & Bundy, K. 2005, 

ApJ, 635, 1006 
Shapley, A. E., Staidal, C. C, Adalbergar, K. L., at al. 

2001, ApJ, 562, 95 
Shen, S., Mo, H. J., White, S. D. M., et al. 2003, MNRAS, 

343, 978 

Somerville, R. S., Hopkins, P. F., Cox, T. J., Robertson, 
B. E., & Harnquist, L. 2008, MNRAS, 391, 481 

Stark, D. P., Swinbank, A. M., Ellis, R. S., et al. 2008, 
Nature, 455, 775 

Steidcl, C. C, Adclbergar, K. L., Shapley, A. E., at al. 

2003, ApJ, 592, 728 

Stewart, K. R., Bullock, J. S., Wachsler, R. H., Mailer, 
A. H., & Zentnor, A. R. 2008, ApJ, 683, 597 

Tacconi, L. J., Ncri, R., Chapman, S. C, et al. 2006, ApJ, 
640, 228 

Teplitz, H. I., McLean, I. S., Becklin, E. E., et al. 2000, 

ApJL, 533, L65 
Tornatore, L., Borgani, S., Dolag, K., & Matteucci, F. 2007, 

MNRAS, 382, 1050 
Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 

2004, ApJ, 613, 898 

Weiner, B. J., Coil, A. L., Prochaska, J. X., et al. 2008, 
ArXiv c-prints 

Weiner, B. J., Papovich, C, Bundy, K., et al. 2007, ApJL, 

660, L39 

Warner, M. W., Roallig, T. L., Low, F. J., et al. 2004, ApJS, 
154, 1 

Wolf, C, Bell, E. F., Mcintosh, D. H., at al. 2005, ApJ, 
630, 771 

Wooslay, S. E. & Weaver, T. A. 1995, ApJS, 101, 181 
Zheng, X. Z., Ball, E. F., Papovich, C, et al. 2007, ApJL, 

661, L41 



